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LONG-TERM  GOALS 

The  long-tenn  goal  of  this  effort  is  to  advance  our  quantitative  understanding  of  the  factors  affecting 
signal  propagation  in  the  marine  environment,  essential  for  radio  tracking,  communication  and 
guidance  applications.  A  significant  issue  in  this  scientific  area  is  the  reproducible  tendency  of  models 
for  propagation  to  overestimate  the  signal's  intensity  at  the  receiver  (Barrios  and  Patterson  (2002)). 
Currently  employed  algorithms  rely  on  the  return  signal's  intensity  for  determining  the  distance  to  an 
object,  thus  a  signal  misinterpretation  has  a  potentially  far-going  practical  consequences.  The  situation 
suggests  that  physical  mechanisms  or  experimental  circumstances  responsible  for  signal  degradation, 
contraction  of  the  coherence  radius,  etc.,  are  not  fully  understood  and  accounted  for. 

A  very  important  phenomenon  influencing  radio  signals  is  the  formation  of  evaporation  ducts 
(refractive  waveguides)  in  the  first  tens  of  meters  over  the  sea  surface.  Models  for  propagation  over  the 
ocean  that  are  currently  in  use,  rely  on  an  averaged  mesoscale  description  of  the  marine  atmospheric 
boundary  layer.  The  only  possible  answer  to  the  intensity  discrepancy  that  could  be  sought  within  this 
approach  is  the  possiblity  of  a  “leaking”  evaporation  duct,  the  possibility  that  the  wavegude  lets 
radiation  escape  into  space.  Two  other  factors,  namely  the  fluctuating  refractivity  in  the  boundary 
layer  and  the  scattering  by  the  ocean  surface  are  generally  ignored.  An  analysis  conducted  within  this 
effort  detemined  that  an  error  in  a  widely  adopted  model  for  scattering  (Miller  et  ah,  (1984))  is, 
enetirely  or  parially,  responsible  for  the  discrepancy  between  model  predictions  and  observational 
results.  Also,  the  complex  motion  of  the  atmosphere  can  be  responsible  for  broadening  of  the 
propagating  beam,  thus  reducing  its  intensity.  A  recently  collected  data  have  produced  clear  evidence 
that  the  boundary  layer  motion  consists  of  both  turbulence  and  wave  effects.  Although  the  role  of  the 
turbulence  has  been  studied  extensively,  the  role  of  the  wave  effects  is  largely  unknown. 

In  this  context,  a  goal  of  this  work  is  to  explore  the  influence  on  signals  of  factors  and  processes  in  the 
marine  boundary  layer  that  have  so  far  been  ignored  (e.g.  micrometeorological  fields).  A  special 
attention  is  given  to  the  multifaceted  role  of  surface  waves,  one  distinct  element  of  the  marine 
environment. 

OBJECTIVES 

The  specific  objective  of  this  effort  is  to  quantify  the  effect  of  the  surface  waves,  manifested  through  a 
modulation  of  the  atmospheric  motion  and  through  defining  the  statistics  of  the  sea  surface,  on  the 
propagation  pattern  of  a  signal  in  the  marine  atmospheric  boundary  layer.  We  seek  a  repersenataion  of 
our  findings  in  a  form  suitable  to  incorporate  into  numerical  models  for  signal  transmission.  Invoking 
the  profound  physical  similarities  between  the  electromagnetic  waves  over  the  ocean  and  the  sound 
signals  in  the  water,  we  also  consider  extending  the  application  of  our  results  to  the  acoustic  case.  The 
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ultimate  objective  is  to  eliminate  particular  contributor  to  the  discrepancy  between  measurements  and 
modeling  results  and  thus  make  numerical  models  for  propagation  more  accurate. 


CSAT 
LI  7500 
Met3A 


CSAT 
LI  7500 
Met3A 


CSAT 

LI7500 


Figure  1.  The  Air-sea  interaction  tower  with  the  instruments  array  during  the  CBLAST  experiment. 
Measurements  of  wind  velocity,  atmospheric  temperature  and  humidity,  pressure  fluctuations,  and 
surface  elevation  were  conducted,  among  others.  Photograph  courtesy  of  Dr.  Jim  Edson. 


APPROACH 

The  data  collected  during  the  CBLAST  field  experiment  (Figure  1)  exhibit  extensive  intervals  (about 
20  %  of  the  experiment's  duration)  of  low  winds.  This  atmospheric  flow  regime  in  these  conditions  is 
illustrated  in  Figure  2.  At  these  low  wind  speeds  the  shear-driven  turbulence  has  low  intensity  as  well 
as  the  turbulent  mixing.  Such  circumstances  allow  vertical  gradients  of  atmospheric  humidity  (and 
refractivity)  to  be  formed  and  maintained. 

Figure  2  demonstrates  that  in  low  wind  conditions  the  motion  of  the  whole  atmospheric  boundary  layer 
is  dominated  by  the  wave  driving.  The  essence  of  mechanism  to  explain  such  wave  modulation  is  the 
fact  that  the  surface  motion  displaces  column  of  air  above  it  and  the  mean  flow  streamlines.  If  the  air 
column  has  a  stratified  distribution  of  the  humidity  (refractivity),  i.e.  in  the  case  of  vertical  humidity 
gradient,  an  instrument  at  a  fixed  height  will  register  humidity  (refractivity)  fluctuations  resulting  from 
the  displacement.  We  will  assume  that  the  surface  waves  have  a  small  slope,  so  the  wave-induced 
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fields  are  linear  responses  to  the  surface  motion  (Hristov  et  al.  (2003)).  This  physical  picture  will  be 
tha  basis  of  our  approach  to  describing  the  structure  of  the  boundary  layer  motion  and  the  influence  of 
the  wave  modulation  on  the  signal  propagation  pattern. 


Figure  2.  At  low  wind  speeds  the  wave  driving  of  the  atmospheric  boundary  layer  is  emphasized. 
In  the  left-hand  column,  the  plots  show  100s  of  the  measured  (not  processed)  signals  of  the  along- 
wind  velocity  at  4  levels  from  the  surface  (top  plot,  colors  indicating  the  height  of  the  instrument 
from  the  surface  in  the  order  blue  -lowest,  green,  red,  cyan  -  highest),  vertical  wind  velocity  (second 
plot,  same  color-height  correspondence),  atmospheric  pressure  fluctuations  (third  plot,  green-lower, 
blue-higher  instrument),  and  surface  elevation  (bottom  plot).  The  plots  in  the  middle  column  show 
the  coherence  between  the  atmospheric  refractivity  and  the  waves  (top  plot,  blue-lowest,  green,  red- 
highest  instrument),  and  the  surface  wave  spectrum  (bottom  plot).  On  the  right-hand  side, 
propagation  regimes  through  idealized  periodic  refractive  media. 


WORK  COMPLETED 

The  propagation  in  an  idealized  periodic  refeactive  media  can  be  illustrated  by  solving  the  ray-tracing 
equation.  The  solutions  (Figure  2)  show  that  the  beam  can  be  confined  along  the  crests  of  the 
refractivity  or  accumulate  displacement  when  propagating  in  an  oblique  direction,  thus  forming 
mirages.  However,  considering  that  the  surface  waves  are  not  monochroomatic,  solving  realistic 
problems  of  propagation  requires  a  statistical  description  of  the  wave-induced  structure  of  the 
boundary  layer.  Important  characteristics  of  the  recieved  signals,  such  as  statistics  of  the  phase 
fluctuations,  of  the  intensity  or  of  the  angle  of  arrival  depend  on  characteristics  of  the  media,  such  as 
the  characteristic  function  of  the  velocity  field,  characteristic  function  of  the  two-point  correlations  of 
the  velocity  field,  structure  function  of  the  atmospheric  refractivity  (Ishimaru  (1978);  Tatarskii 
(1971)).  We  evaluated  these  functions,  as  follows. 
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We  assume  that  the  velocity  field  in  the  atmospheric  boundary  layer  u  can  be  decomposed  into  mean 
u  ,  turbulent  u  and  wave  induced  u  components  as  u  =u  +ii'  +  u  .  Within  the  asumption  of  small- 
slope  waves  ( ka  « 1 )  the  wave-induced  field  u  =  {/7,v}  ( u  being  the  along- wind  component  of  the 
wave-induced  velocity  and  v  being  the  vertical  component)  can  be  considered  as  a  linear  response  of 
the  surface  wave  forcing  77  =  ae~lk(x~ct) ,  i.e.  {u,  v  }  =  \tf)u ,  tj)v }//  =  (u„  /  k){-  (dtp  /  dy),  ik(p\rj ,  where  tp  is  a 
solution  of  the  Rayleigh  equation  (Hristov  et  al.  (2003)). 

Characteristic  function  of  the  atmospheric  velocity.  The  characteristic  function  of  the  atmospheric 
velocity  /(k  r)  =  (exp(ik  ■  vr)j ,  through  the  radar  equation  (Ishimaru  (1978))  detennines  the 

correlation  function  of  the  signal’s  electrical  field  as  well  as  the  correlation  function  of  the  signal’s 
phase  (Wheelon  (2001)).  Our  purpose  here  is  to  incorporate  the  wave  effects.  Invoking  the 

decomposition  u  =  u  +  u'  +  u  and  considering  the  turbulence  u  and  the  wave  effects  u  =  {u, v }  to  be 
statistically  independent,  we  have  shown  that  the  characteristic  function  factorizes  into  a  turbulent 
Zrfkr)  and  wave-induced  /fk  r)  multipliers: 

Zu(kT)  =  Zu'(kT)Zz(kT) . 


To  model  the  wave-induced  part, we  employ  two  considerations.  First,  that  the  ocean  surface  can  be 
viewed  as  an  ergodic  random  process  which  is  a  superposition  of  independent  Fourier  modes  with 
random  phases.  Relying  on  the  central  limit  theorem,  we  can  reasonably  conjecture  that  the  ocean 
surface  has  a  Gaussian  distribution,  both  in  space  and  time.  Second,  recalling  our  assumption  that  the 
wave-induced  fields  {u,v}  =  {(f ,  (f  }>]  =  (u*  /  k){- (dtp  I  dy),  iktp\rj  are  linear  responses  to  the  wave 
forcing,  and  invoking  the  theorem  that  a  linear  transform  preserves  that  Gaussianity  of  a  random 
process  (Parzen  (1962)),  we  can  conclude  that  the  wave-induced  velocities  u  =  {«7,  v}  also  should  be 
Gaussian  random  variables.  Consequently, 


p(u)  =  (2^cr|)”3/2  exp 


and  the  characteristic  function  for  the  wave-induced  velocity  is  then 

Xu  (k  t)  =  exp  [-  {k2trl  r2 )/  2] . 


The  standard  deviation  is  obtained  considering  ut  =  <7~  +  a~  and  that 


=XX\dc=(ujK)xX}rdc 


where  u*  is  the  friction  velocity  and  S  =  (////  ')  is  the  surface  wave  spectrum 
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The  characteristic  function  of  the  two-point  difference  of  the  atmospheric  velocity 


Za  «w2)(£0  =  (exp[/£r(wfoO-K(r2,f))J_ 

is  needed  to  calculate  the  correlation  function  of  the  signal’s  intensity  fluctuations  (Tatarskii  (1971)). 
For  Gaussian  processes  JAl7(j5)(kr)  =  exp 


-  j ikzf  ({(k  /  k )  •  A u ( p ,  t)J 


with 


(( k  /  k )  •  Au(p,  t)J  \  =  -j kfjDy  ip)  and  Dy  (p,  r,  t )  =  (jul  (r  +  p,t )  -  ui  ( r ,  t)\uj  (. r  +  p,t )  -  Uj  (r ,t)]) . 


Again,  considering  the  decomposition  u  =  u  +  u  +  u  and  that  the  turbulence  u  and  the  wave  effects 
u  =  {/7,v}  are  uncorrelated  and  statistically  independent,  the  structure  function  splits  into  turbulent 
D'j(p,r,t)  and  wave-induced  D:/(p,r,t)  parts  as 

Dtj  (A  r,t)  =  D'^  (p,  r,t)  +  Dy  (p,  r,  t ) 


and  the  characteristic  function  factorizes  into  turbulent  and  wave-induced  multiplyers 
XmpfiT)  =  [%mP)ikr)\x/,irp)(kT)_ . 


It  is  important  to  note  that  while  for  isotropic  turbulence  the  turbulent  structure  function  (as  well  as  the 
corresponding  characteristic  function)  depends  only  on  p  =  ^  -  r2 1 ,  their  wave-induced  counterparts 
are  anisotropic.  Here  again,  we  can  use  the  solutions  of  the  Rayleigh  equation  (Hristov  et  al.  (2003)) 
{u,  v }  =  ,  (j)~  }rj  =  (w*  /  at){-  ( d(p  /  dy),  ik(p}r/  to  find  the  explicit  fonn  of  r  3  (k  r) . 

AU(p) 


The  structure  function  of  the  atmospheric  refractivity  Dn  (q ,  r2)  =  (jn(i[ )  -  n(r2  )]2  ^  occurs  when  one 

needs  to  quantify  the  statistics  of  the  fluctuations  of  signal’s  phase  and  angle  of  arrival.  Assuming  that 
the  fluctuations  of  atmospheric  refractivity  over  the  ocean  n  consist  of  turbulent  n  and  wave-induced 
n  parts  n  =  n'  +  n  and  that  n  and  n  are  uncorrelated,  i.e.  (n'n'j  =  0  ,  the  structure  function  splits  into 

turbulent  D’n(T\,r2)  ==  ^n'(fx) -n'(r2)]2^  and  wave-induced  parts  Dn(i\,r2)  =  ^[«(q)-n(F2)]2^  as 


Dn  (f  A)  =  D'n  (h  A)  +  Dn(ff2). 


The  turbulent  structure  function  D'n(T\,r2)  has  been  studied  extensively  since  the  1960s  (Tatrskii 
(1971)),  yet  nothing  is  known  about  the  wave-induced  part  Dn  if ,  r2 ) .  Unlike  D'n  (f  ,f)  =  D'n  (|q  -  f  |) , 

the  wave  induced  tenn  Dn(f,r2)  of  the  structure  function  is  anisotropic,  since  the  vertical  direction 
plays  a  distinctly  different  role  than  the  two  horizontal  directions. 
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To  evaluate  Dn(rl,r2)  we  will  assume  that  the  waves  are  displacing  vertically  the  mean  flow 

streamlines  and  the  column  of  air  with  stratified  refractivity.  The  change  in  refractivity  registered  by 

an  instrument  at  a  fixed  height  is  expressed  through  the  vertical  gradient  of  the  refractivity 

(dN (z)  /  dz)  and  the  vertical  displacement  of  the  streamline  at  height  Z  caused  by  a  spectral  mode  in 

rdN(Z)^ 


the  wave  spectrum  with  phase  speed  c  ,  &  =  J  v  ( Z,k)dt 


as  n  =  ■ 


dZ 


(dz)  =  T (Z,k)r/(k) .  We  have 


introduced  T(Z.k)  as  a  transfer  function,  relating  the  surface  displacement  //(c)  and  the  refractivity 
fluctuations  n  .  From  here,  introducing  r  =  R  +  Ze_  =  Xex  +  Yev  +  Ze _ ,  the  correlation  function  of  the 
refractivity  can  be  expressed  as 


d(Rx,Zx)d(R2,Z2))  =  J  /kM)T(zJ)T\Z2,k)Sm(k)dk 


which  is  sufficient  to  determine  also  (n(R,Z)n(R,Z)j  and  therefore  the  wave-induced  tenn  of  the 
structure  function  Dn (T\ ,  r2 ) . 

RESULTS 

Numerical  models  for  propagation  in  ducting  conditions  have  shown  a  persistent  tendency  to 
overestimate  the  intensity  at  the  receiver  (Barrios  and  Patterson  (2002)),  thus  leading  to  overestimation 
of  the  distance  to  the  object  producing  the  radar  return  as  well  as  to  overestimation  of  the  available 
response  time.  We  reviewed  the  possible  mechanisms  for  signal  degradation,  potentially  responsible 
for  this  deficiency  in  propagation  models.  We  determined  that  the  widely  adopted  model  for  sea 
surface  scattering  (Miller  et  al.  (1984))  contains  an  error,  which,  entirely  or  partially,  leads  to  the 
observed  error  in  propagation  modeling.  We  proposed  and  alternative  of  the  model  of  Miller  et  al. 
(1984),  which  is  free  of  that  model  deficiencies.  We  explored  the  possibility  that  the  atmospheric 
motion  is  responsible  for  broadening  of  the  radar  beam,  thus  effectively  reducing  its  intensity  and 
leading  to  signal  degradation.  Data  from  the  CBLAST  experiment  have  indicated  that  in  certain 
conditions  wave  modulation  can  dominate  the  motion  in  the  atmospheric  boundary  layer.  We 
quantified  the  influence  of  that  modulation  on  the  propagation  pattern  by  expressing  the  characteristic 
function  of  the  atmospheric  velocity,  the  characteristic  function  of  two-point  differences  and  the 
structure  function  of  the  atmospheric  refractivity  in  a  form  suitable  for  incorporating  in  numerical 
propagation  models. 

IMPACT/APPLICATIONS 

The  analysis  outlined  above  has  reviewed  and  identified  several  causes  for  discrepancy  between 
propagation  model  results  and  observations.  The  proposed  corrections,  in  describing  the  scattering  by 
the  rough  ocean  surface  and  by  quantifying  the  influence  of  the  atmospheric  boundary  layer  motion  on 
the  propagation  pattern,  are  expected  to  substantially  improve  the  performance  of  models  for  signal 
propagation  over  the  ocean.  Although  the  strong  modulation  of  atmospheric  motion  has  been  observed 
over  the  mid-latitude  Atlantic,  it  is  likely  that  the  phenomenon  would  be  both  more  prevalent  and  more 
pronounced  over  the  Indian  Ocean,  where  low  wind  conditions  are  often  encountered  and  where  dry 
air  from  the  surrounding  deserts  can  move  over  the  ocean  and  cause  strong  refractive  ducts.  Because  of 
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the  profound  physical  similarities  between  atmospheric  electromagnetics  and  underwater  acoustics,  the 

application  of  these  results  can  be  extended  to  sound  signals  in  the  ocean. 

RELATED  PROJECTS 

The  PI  is  unaware  of  any  related  ONR  sponsored  projects. 
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